Ground-state phase diagram of a spin-| frustrated ferromagnetic XXZ chain: 
Haldane dimer phase and gapped/gapless chiral phases 
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The ground-state phase diagram of a spin-^ XXZ chain with competing ferromagnetic nearest- 
neighbor (Ji < 0) and antiferromagnetic second- neighbor (J2 > 0) exchange couplings is studied 
by means of the infinite time evolving block decimation algorithm and effective field theories. For 
the SU(2)-symmetric (Heisenberg) case, we show that the nonmagnetic phase in the range —4 < 
J1/J2 < has a small but finite ferromagnetic dimer order. We argue that this spontaneous dimer 
order is associated with effective spin-1 degrees of freedom on dimerized bonds, which collectively 
form a valence bond solid state as in the spin-1 antiferromagnetic Heisenberg chain (the Haldane 
spin chain). We thus call this phase the Haldane dimer phase. With easy-plane anisotropy, the 
model exhibits a variety of phases including the vector chiral phase with gapless excitations and the 
even-parity dimer and Neel phases with gapped excitations, in addition to the Haldane dimer phase. 
Furthermore, we show the existence of gapped phases with coexisting orders in narrow regions that 
intervene between the gapless chiral phase and any one of Haldane dimer, even-parity dimer, and 
Neel phases. Possible implications for quasi-one-dimensional edge-sharing cuprates are discussed. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.80.+q 



I. INTRODUCTION 

The search for novel quantum states in frustrated mag- 
nets has been a subject of intensive theoretical and ex- 
perimental research. Onc-dimcnsional (ID) systems offer 
unique laboratories for this search, as strong fluctuations 
enhance the tendency toward unconventional quantum 
states.El Among them, the ID XXZ model with competing 
nearest-neighbor J\ and second-neighbor Ji interactions, 
defined by the Hamiltonian 



(1) 



provides a paradigmatic example expected to host rich 
variety of physics. Here Si = (Sf, S^,Sg) represents the 
spin-i operator at the site I £ Z and A parametrizes the 
XXZ exchange anisotropy. The model has frustration as 
far as J2 is antiferromagnetic, irrespective of the sign of 
Ji. 

Early theoretical studies on the model ([j]) mostly 
considered the pA§e when both Ji and J2 are 
antiferromagneticialD However, interest is now growing in 
the case of ferromagnetic J\ < and antiferromagnetic 
J2 > because of its relevance to quasi-lD edge-sharing 
cuprates. Among such cuprates, LiCu 2 02 (Rcfs. BO), 
LiCuV0 4 (Rcfs. [fijjn]), and PbCuS0 4 (OHk^cf.@, 
for example, exhibit multifcrroic behaviorSjtSo i.e., spi- 
ral magnetic orders and concomitant ferroelectric polar- 
ization at low temperatures. The negative sign of J\ 
indeed plays a key role in stabilizing tha vector chiral 
order responsible for these phenomcna.E-3 By contrast, 
Rb 2 Cu2Mo30i2 (Ref. |l6|) shows no sign of magnetic or- 
der down to very low temperatures and may be consid- 



ered as a candidate system for a spin liquid or a valence 
bond solid. 

In this paper, we study the ground-state properties of 
the spin-i frustrated ferromagnetic XXZ chain (Q) with 
J 1 < and J 2 > 0, by means of the infinite time evolving 
block decimation algorithm (iTEBD)E-D and effective field 
theories based on the bosonization methods. Previous 
works on the case with easy-plane anisotropy < A < 1 
have discussed the competition among the vector chiral 
phase with gapless excitations ajid Ihe^dimer and Neel 
phases with gapped excitations JjH'EtEj The main goal 
of this paper is to present a conclusive phase diagram of 
the model ([j]) , which is shown in Fig. [j], thrausjh-dctaired 
analyses that extends our previous works .EjoEj Firstly, 
we uncover the nature of the nonmagnetic phase around 
the SU(2)-symmetric case A = 1, which has long been 
controversial. We show that this phase has a dimer order 
associated with an emergent spin-1 degree of freedom on 
every other bond. We term this new phase the Haldane 
dimer phase. Secondly, we show the existence of narrow 
gapped phases that intervene between the gapless chiral 
phase and any one of gapped dimer and Neel phases. As 
weak inter-chain couplings are turned on, while the gap- 
less chiral phase evolves into a spiral magnetic order, the 
Haldane dimer phase can be stabilized by a coupling with 
phonons due to the spin-Peierls mechanism. Our phase 
diagram may thus provide a useful starting point for un- 
derstanding the competing phases in quasi-lD cuprates. 

Let us briefly review previous results on the model (jl|) 
and summarize our new findings. While we are mainly 
concerned with the case of J\ < and J2 > in this 
paper, for comparison, we also review established results 
on the case of antiferromagnetic J\ , J2 > alongside. 

In the classical limit S — > 00, the ground state phase 
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FIG. 1: (Color online) The ground-state phase diagram of 
the model (|l]) with Ji < and J2 > 0. The vector chiral 
phase, which has a non- vanishing vector chirality (H) , extends 
between the two boundaries with the "□" symbols. These 
boundaries are determined as in Figs. || and |5| (see the verti- 
cal solid lines in these figures). Around the highly degenerate 
point (Ji/ J2, A) = (—4,1), these two boundaries could not 
be determined accurately, but we expect both of them to con- 
tinue to this point. It was also difficult to draw the boundary 
around the right top corner of the phase diagram. The onsets 
of (Haldane and even-parity) dimer and Neel orders occur in- 
side the vector chiral phase, as indicated by the "x" symbols 
(determined as in Figs. and [|). Thus there are narrow in- 
termediate phases (between the "□" and "x" symbols) where 
two kinds of orders coexist. The phase boundaries among the 
TLL, even-parj-fev dimer, and Neel phases are determined in a 
previous work.Ej On the right of the "O" symbols, the even- 
parity dimer and Neel phases alternately appear when ap- 
proaching the point (J1/J2, A) = (—4, 1): the first transition 
occurs at A ss-iO.7 ("o" symbols) and the second at A w 0.93 
(not shown) .cj On the line where A = 1 and J1/J2 < —4, 
the ground state is ferromagnetic (FM). The "+" symbols in- 
side the dimer and Neel phases indicate the Lifshitz line, on 
which the short-range spin correlation changes its character 
from incommensurate (IC) to commensurate (C); see Fig. 



diagram of Eq. ([!]) does not depend on A in the range 
< A < 1. The ground state has ferromagnetic or- 
der for J1/J2 < —4 and antifcrromagnctic (Neel) or- 
der for J1/J2 > 4. For < \J\\/ J2 < 4, the ground 
state is in a spiral magnetic ordered phase, in which 
the spins rotate by an incommensurate pitch angle Q = 
± arccos(— J1/4J2) along the spin chain. Except for the 
isotropic case A = 1, the spiral plane is fixed in the xy 
plane, and the vector chirality 



^1,1+1 
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has a non- vanishing uniform value nf e+1 = ± sinQ inde- 
pendent of I. Here (■ ■ • ) stands for average in the ground 



state (with long-range order, if any). 

In the ground state of the quantum spin-i model, a 
long-range magnetic order with broken U(l) spin rota- 
tional symmetry is generally prohibited, unless the uni- 
form magnetic susceptibility is divergent as in the case 
of ferromagnetism.L3 However, a long-range order (LRO) 
of the vector chirality k\ l+1 that breaks only the Z2 
parity symmetry can survive quantum fluctuations in 
the case of A 7^ 1. Using the bosonization theory for 
\Ji\/J2 *C 1 and < A < 1, Nersesyan et alU pre- 
dicted the appearance of the vector chi ral pha se with 
gapless excitations (as reviewed in Sec. IVA2). This 
gapless chiral phase shows the spatially uniform vector 
chirality k\ e+1 7^ and power-law decaying (incommen- 
surate) spiral spin correlations; this phase may therefore 
be viewed as a quantum counterpart of the classical spi- 
ral phase. The gapless chiral phase competes with other 
quantum phases, in particular, valence bond solids driven 
by quantum fluctuations. In fact, for antifcrromagnctic 
Ji > 0, a dimerized phase, in which the singlet state 
(I'M) - I 4-T ) )/v / 2 (written in the {Sf} basis) is formed 
on dimerized bonds, appears in -a .large part of the classi- 
cal spiral regime < J\j J 2 < 4j3l3 and the gapless chiral 
phase appears only in a small region in the space spanned 
by J 1 /J 2 and A.Q 

The phase diagram for the case of ferromagnetic J\ < 
and easy-plane anisotropy < A < 1 is presented in 
Fig. |l]. Early workalao mainly discussed the transition 
from the Tomonaga-Luttinger liquid .XTLL) phase to a 
dimer phase with an even-parity unitcj | f-i ) +4 i-t ) ap- 
pearing for < A < 0.7. Our recent workal3'Ej have 
uncovered a rich phase structure in an extended param- 
eter space of J1/J2 and A. In Rcf. 15, it was shown that 
the gapless chiral phase appears in a wide region, and sur- 
vives up to the close vicinity of the isotropic case A = 1 
for -4 < Ji/ J 2 < -2.5 (we also refer to Refs. |||| for 
related earlier works). This remarkable stability of the 
gapless chiral phase for J\ < indicates that the sign of 
J\ plays a crucial role in stabilizing the vector chirality 
and the associated ferroelectric polarization in multifcr- 
roic cupratcsBrJ In Rcf. [24|, the instability of the TLL 
phase toward gapped phases was analyzed using the ef- 
fective sine-Gordon theory combined with numerical di- 
agonalxzation. It was found that the even-parity dimer 
phaseEJ discussed in Refs. [l8|-p0|and a Neel ordered phase 
appear alternately as A is increased on the right side of 
the TLL phase in the phase diagram. 

An important result of this paper is concerned with 
the nature of the nonmagnetic phase for —4 < J1/J2 < 
around the SU(2)-symnic±ric case A = 1. Previous 
field-theoretical analyseaJ'EJ have suggested that a dimer 
phase with a very small energ y gap should appear in 
this region (as reviewed in Sec. [II A). However, neither 
a dimer order nor an energy gap has been detected in pre- 
vious numerical studies. Using the iTEBD, which allow 
us to treat infinite-size systems directly, we present the 
first numerical evidence of a finite dimer order parameter 
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Remarkably, this dimcr order is associated with ferro- 
magnetic correlations (Se ■ SV+i) > of alternating 
strengths, in contrast to antifcrromagnctic correlations 
in singlet dimers for J\ > 0. In this case, it is natu- 
ral to interpret that effective spin-1 degrees of freedom 
emerge on the bonds with stronger ferromagnetic corre- 
lation, forming a valence bond solid stateEa as in the Hal- 
dane chainO We thus call this new phase the Haldanc 
dimer phase. 

We also present detailed analyses of the_aiiisotropic 
case A ^ 1, extending our previous works. Ej'lII In par- 
ticular, we analyze the transition from the gapless chi- 
ral phase to each of the Haldane dimer, even-parity 
dimer, and Neel phases, and identify narrow interme- 
diate gapped phases where two kinds of orders coexist 
(the regions between "□" and "x" symbols in Fig. [I]). 
Furthermore, we describe how the properties of various 
phases can be,|Gaptured in the language of the Abelian 
bosonizationMEll for \J\\/ J2 "C 1, as summarized in Ta- 
ble I 

The rest of the paper is organized as follows. In Sec. ||, 
we present the numerical results on the order parameters 
and half-chain entanglement entropy, which provide the 
most basic information for identifying symmetry-broken 
phases. In Sec. Ill, we discuss in detail the dimer phases 
in the SUl2)pSymmetric case A = 1 from both field- 
theoretical 



and numerical analyses. In Sec. IV 



discuss the case with the easy-plane anisotropy < A < 
1. In particular, we review the effective field theory for 
the gapless chiral phase, U and, following Ref. |34|, discuss 
its instability towards gapped chiral phases. The ranges 
of the gapped chiral phases are then determined numer- 
ically by analyzing the spin correlations. In Sec. |v|, we 
brie fly , des cribe how the quantum phases in the easy-axis 
caseE3E3c3 can be understood in the Abelian bosoniza- 



tion framework. In Sec. VI, we conclude the paper and 
discuss implications of our results for quasi-lD cuprates. 



II. NUMERICAL ANALYSIS OF ORDER 
PARAMETERS 

In this section, we present numerical results on several 
order parameters and half-chain entanglement entropy 
calculated by iTEBD. The vector chiral order parame- 
ter and the entanglement entropy are used to determine 
the boundaries of the region where the long-range vec- 
tor chiral order exists (the "□" symbols in Fig. [j]). The 
numerical results in this section also suggest the exis- 
tence of the narrow intermediate phases (between "□" 
and "x" symbols) in which the vector chiral order coex- 
ists with the dimer or Neel order. The precise ranges of 
these intermediate phases, however, will be determined 



in Sec. IV B 2 



Before presenting the numerical results, let us briefly 
note characteristic features of our numerical method; for 
more detailed account of the method, see Supplementary 
Material of Ref. pi. The iTEBD algorithmic we em- 



ployed is based on the periodic matrix product represen- 
tation of many-body wave functions of an infinite system. 
It can directly address physical quantities in the thermo- 
dynamic limit, and is free from finite-size or boundary 
effects. The (variational) wave function is optimized to 
minimize the energy. The precision of the algorithm is 
controlled by the Schmidt rank x, which gives the linear 
dimension of the matrices. We exploited the conservation 
of the total magnetization J^e Sf = to achieve higher 
efficiency and precision of the calculations. When this 
algorithm is used in ordered phases, a variational state 
finally converges to a symmetry-broken state with an as- 
sociated finite order parameter (if it is.-allowed by the 
periodicity of the matrix product state) S3 In our imple- 
mentation, we used a period-4 structure for the varia- 
tional matrix product state. In this setting, the vector 
chiral, dimer, and Neel order parameters analyzed in this 
section can all be calculated through local quantities. In 
order to allow a finite vector chiral order parameter, the 
initial state must contain complex elements "seed" 
for the symmetry breaking £3 



A. Vector chiral order 

Figures || and ^ present our numerical results along 
the vertical line 3\j Ji = —2 and the horizontal line 
A = 0.8, respectively, in the phase diagram (Fig. [j]). 
Let us first look at the vector chiral order parameter 
n\ 2 = ((Si x S 2 ) z ) displayed in Figs. |(a) and |(a). This 
order parameter is always found to be spatially uniform 
along the spin chain in the present model, so we have 
fixed the site labels. By observing the rapid increase of 
Ki 2 , we find the onset of the vector chiral phase. It is 
natural to think that this rapid increase comes from the 
Ising nature of the transition with exponent /? = 1/8 for 
the spontaneous order parameter, as previously demon- 
strated in the XY case A = 0.E3 To determine the transi- 
tion points more precisely, however, we use the half-chain 
entanglement entropy explained next. 

The half-chain van Neumann (vN) entanglement en- 
tropy is defined asliZI 



5 V1N = - J2 A " 

a=l 



In Xz 



(4) 



where {X a } is a set of Schmidt coefficients associated 
with the decomposition of the infinite system into the 
left and right halves and x 1S t ne Schmidt rank. As the 
system approaches a critical point characterized by a con- 
formal field theory witb,£ucentral charge c, this quantity 
is known to diverge asL3c9 



G 



ln£ + si, 



(5) 



where £ is the correlation length and Si is a non- universal 
constant. In an iTEBD calculation with a finite Schmidt 
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FIG. 3: (Color online) (a) Chiral and (b) Neel order parame- 
ters and (c) half-chain entanglement entropy as a function of 
J1/J2 for fixed A = 0.8, calculated by the iTEBD. The solid 
and broken vertical lines indicate the onsets of the vector chi- 
ral and Neel orders, respectively. The former is determined 
by the peak position in the entanglement entropy in panel 
(c), while the latter is determined in Fig. |l5|. In the narrow 
intermediate phase in —3.225 < J1/J2 < —3.200, the vector 
chiral and Neel orders coexist. 



FIG. 2: (Color online) (a) Chiral and (b,c) dimer order param- 
eters and (d) half-chain entanglement entropy as functions of 
A for fixed J1/J2 = —2. These are calculated by the iTEBD 
with Schmidt ranks \ = 200 and 300. Panel (c) is a zoom 
of panel (b). In panel (d), Si™ an d £2,3 are defined for the 
bipartitions of the system at the bonds (1,2) and (2,3), re- 
spectively. Solid vertical lines indicate the boundaries of the 
vector chiral phase, and are determined from the onsets of 
the vector chiral order parameter in panel (a) or more accu- 
rately from the peaks in the entanglement entropy in panel 
(d). Broken vertical lines indicate the transition points on 
which dimer orders set in. These points are difficult to locate 
within the analysis of dimer order parameters in panels (b,c), 
and are instead determined by the analysis of spin correlations 
in Fig. [l4|. Narrow intermediate phases exist between solid 
and broken vertical lines, where the vector chiral and dimer 
orders coexist. In the intermediate phase in 0.61 < A < 0.63, 
a dip in the entanglement entropy is seen in panel (d), as 
zoomed in the inset. 



rank x, the divergence of S vN at the critical point is 
replaced by the increasing function of Xf-^ 



5 vN = 



y/U/c+1 



(6) 



where is another non- universal constant. The calcu- 
lated entanglement entropy is shown in Figs. §(d) and 
"|(c). In Fig. ||(d), we plot two entropies and SjSJ as- 
sociated with the bipartitions of the system at the bonds 
(1,2) and (2, 3), since these bonds are inequivalent in the 
neighboring dimer phases. By finding peaks of S' vN , we 
can determine the boundaries of the vector chiral phase, 
more accurately than by using k| 2 ; see the solid vertical 
lines in Figs. ^| and [| In this way, we have determined 
the square symbols in Fig. [j]. Although we could not ex- 
tract c from the current data of 5* vN using Eq. (0) (which 
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is expected to be satisfied for larger x), it is natural to 
expect that these critical points are characterized by the 
two-dimensional Ising universality class with c = 1/2 (we 
again note that the critical exponent /3 = 1/8 for this 
class was confirmed in the XY caseEa). 

In most part of the vector chiral phase, the entangle- 
ment entropy increases as a function of %, indicating a 
critical nature. .-Indeed, in the effective field theory of 
Nersesyan et a/., El the gaplcss chiral phase has c = 1, and 
the increase of S from the cases of x = 200 to 300 is 
roughly consistent with AS = 0.224 ln(300/200) = 0.091 
expected from Eq. (^|) for c = 1. Near the boundaries 
(solid vertical lines), the entanglement entropy shows 
dips, whose implications will be discussed later. 



B. Dimer orders 

Next we look at the xy and z components of dimer 
order parameters, 
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The alternation of the sign of D^ v l+1 e+2 or Df^ , 1 g+2 
along the spin chain would indicate some sort of dimer 
ordering. We assign the site labels in such a way that 
^123 < 0- The two order parameters are plotted in 
Figs. g(b,c). We find that -0^23 an d -D123 are both finite 
and have mutually opposite signs for A < 0.65. By con- 
trast, the two order parameters have small finite values of 
the same sign for A > 0.9; in spite of the smallness, they 
are rather stable when the Schmidt rank \ is increased 
as seen in the zoomed plot in Fig. ||(c). These results 
indicate that the dimer phases in the two regions are of 
distinct types. 

The nature of the dimcn-phase for A < 0.6 can be eas- 
ily understood as follows.EJa In the XY limit A = 0, 
the sign of J\ in Eq. ([!]) can be reversed by perform- 
ing the 7r rotations of spins around the z axis on every 
second sites. From the fact that the doubly degenerate 
ground states at (Ji/J 2 ,A) = (2,0) are given by the 
products of singlet dimers, one finds, through the above 
7r-rotation transformation, that the exact ground states 
at (Ji/J2,A) = (—2,0) are given by the dimer states 
whose unit is now replaced by ( | TT ) + I TT ) ) / V% (written 
in the {Sf} basis). We note that this unit has the even 
parity with respect to the inversion about a bond cen- 
ter, in contrast to the odd parity of the singlet dimer at 
J\ > 0. The direct product states of even- parity dimers 
-2Df 23 = ±1/2. The mutually opposite 
and Df 23 an d t ne approximate relation 
D i23 ~ -2£>f 2 3 found for A < 0.6 in Fig. |(b) indicate 
that the even-parity nature of the dimer unit persists in 
this region. We thus call this phase the even-parity dimer 
phase E3 It is distinct from the singlet dimer phase ap- 
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pcaring for J\ > 0, in which D 123 emu -i^ 123 
same sign. 



and Df o show the 



2 3 and D\ 23 are 



In the region A > 0.9 in Fig. g, D*$.. 
both negative as in the singlet dimer phase. However, 
forming nearest-neighbor singlet dimers is unlikely for 
ferromagnetic J\ < 0. In Sec. II], we point out that 



the dimer order in this region is associated with ferro- 
magnetic nearest-neighbor correlations (Se ■ <§^_|_i) > of 
alternating strengths along the chain, in marked contrast 
to an antifcrromagnetic correlation in a singlet dimer. A 
more detailed comparison of the dimer phases for J\ < 
and Ji > in the isotropic case (A = 1) will be presented 
in Sec. EH. 



In the region of a finite vector chiral order (0.61 < 
A < 0.92) in Fig. g, we find that the two dimer order 
parameters remain finite in the narrow regions between 
the solid and broken vertical lines. This indicates the 
existence of the chiral dimer phases (originally predicted 
in Ref. |34|), in which the vector chiral and dimer orders 
coexist and there are four-fold degenerate ground states 
below an excitation gap. In the entanglement entropy, a 
dip is seen in the interval 0.61 < A < 0.63, which also 
supports the existence of an intermediate gapped phase. 
The peaks in the entanglement entropy indicated by the 
solid lines in Fig. g(d) correspond to the Ising critical 
point between two gapped phases. Between the two bro- 
ken lines in Fig. |^, the dimer order parameters diminish 
and the entanglement entropy increases as we increase 
the Schmidt rank these features are consistent with 
the gapless chiral phase. The precise determination of 
the phase boundaries between gapped and gapless chiral 
phases is difficult within the analysis of the order param- 
eters and entanglement entropy in Fig. ||; it will be done 
instead by ana lyzing spin correlation functions in Fig. |IJ 
in Sec. EVB2l ' 



Neel order 



The appearance of a Neel phase with spontaneous stag- 
gered magnetizations {Sf) cx (—1)^ is discussed in detail 
in Ref. |2J. In Fig. ||(b), this Neel order is detected in 
the region J1/J2 < —3.2 by measuring (Si). As in the 
case of the dimer phases, even in the region where the 
vector chiral order is finite (J1/J2 > —3.225), the Neel 
order parameter remains finite. This indicates the exis- 
tence of a narrow chiral Neel phase, in which the vec- 
tor chiral and Neel orders coexist. The ground states in 
this phase should be four-fold degenerate with a finite 
excitation gap. In Fig. ||(c), a dip in the entanglement 
entropy can be found in this region, consistent with the 
expected gapped excitation spectrum. The precise deter- 
mination of the transition point will be done in Fig. 



in Sec. IV B 2 



6 



III. ISOTROPIC CASE A = 1 



In this section, we present detailed analyses of the 
model (|j) in the isotropic case A = 1. While it is 
known that the.|-singlet dimer phase appears for < 
J\j J% ^ 4.15,oHe3 the nature of the nonmagnetic ground 
state in —4 < Ji/J 2 < has not been well understood. 
Sec. " 



In Sec. ilftA, we summarize previous field-theoretical 
analysesoEJfor the weak-coupling limit | J% | <C J 2 , which 
predicted the appearance of dimer orders for both signs 
of Ji . At first glance, this result may seem bizarre since 
the singlet dimerization on the J\ bonds, as formed in 
the case of antiferromagnetic J\ > 0, is unlikely to occur 
in the case of ferromagnetic J\ < 0. In Sec. [II B, we 



present our numerical results and point out a remarkable 
difference between the J\ > and J\ < cases in the 
way how the system hosts the dimer order. This leads us 
to propose the picture of the "Haldane dimer phase" for 
the dimer phase with J\ < 0. Although the ground-state 
wave functions are largely different between the Haldane 
and singlet dimer phases, we argue that the two phases 
in fact share a common hidden order. 



A. Field-theoretical analyses 

Here .-..we . .-summarize previous field-theoretical 
analysesOOcJEJ for |Ji| <C Ji- In this regime, the 
model (0) can be viewed as two antiferromagnetic 
Heisenberg spin chains which are weakly coupled by the 
zigzag interchain coupling J\ as in Fig. [|. We apply 
the Abelian and non-Abelian bosonization techniques to 
describe the two chains separately, and then treat the 
interchain coupling J± as a weak perturbation. 



1. Non-Abelian bosonization 

We start from the non-Abelian bosonizatiorj&El de- 
scription of the isotropic model (|l|) with A = 1, and 
present the renormalization group (RG) analysis to iden- 
tify (marginally) relevant perturbations. 

In the limit J\/J 2 — > 0, each isolated antiferromag- 
netic Heisenberg chain is described by the SU(2)i Wcss- 
Zumino-Witten (WZW) theory, with the spin velocity 
V = (7r/2)J2a, per. tu . rbed by a marginally irrelevant 
backscattering term.EM3c3 The spin operators in the n- 
th chain (n = 1,2) can be decomposed as 



S 2 j+n -> a[M n {x n ) + (-l) j N n (x n )] 



(8) 



with x\(j) = (j — |)a and x 2 (j) = (j + V)a, where a is 



the lattice spacing of each chain; see Fig. |4|. The uniform 
and staggered components, M n and N n , have the scal- 
ing dimensions 1 and 1/2, respectively. The former can 
be decomposed into chiral (right and left) components: 
M n = M n ji + M n L- Another important operator is the 



A ' ; A A — : 
y ■' \ -1 yr \.--„-, 
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FIG. 4: Zigzag chain picture for the J1-J2 chain model (|l|). 
The x axis indicates the coordinate for the continuum descrip- 
tion. 



(in-chain) staggered dimerization operator e„ define by 
{-l) 3 S 2j + n ■ S 2] + n +2 -> ae n (x n ), (9) 

which has the scaling dimension 1/2. 

The inter-chain zigzag coupling J\ produces at most 
marginal perturbations, in the RG sense, around the 
WZW fixed point; relevant perturbations such as Ni-N 2 
are prohibited by the symmetry of the zigzag chain 
model. The symmetry-allowed marginal perturbations 
are summarized as 



H' = J dx^giOi, 



(10) 



where i runs over the following five operators@ 

(11a) 
(lib) 
(11c) 

KJVi.^JVa-JVa-SxJVi), (lid) 



O bs = Mm ■ Mil + M 2R ■ M 2L , 

d = Mifl • M 2L + Mil ■ M 2Rl 

2 = Mm ■ M 2R + Mil ■ M 2L , 
a 



O 



Cdtw = -^(eid x e 2 - e 2 d x ti) 



(lie) 



Here Obs is the backscattering term present in isolated 
chains. The zigzag J\ coupling produces the current- 
current interactions, Oi and 2 , and the twist operator 
O tw . The dimer twist operator Odtw is generated in the 
RG process as we see later. The bare coupling constants 
are given by 

5bs (0) = -0.23(2tto), fli(0)=52(0) = 2Jio, (12) 
flw(0) = Jia, ffdtw(O) = 0, (13) 

where ghJO) was estimated in Ref. |^|. All the operators 
in Eq. (gdj) have the scaling dimensions 2, and their com- 
petition in the RG flow must be analyzed carefully by 
deriving the RG equations. We define the dimensionless 
coupling constants 



Gi = ^ (^bs,l,2), 



Gi = 



2irv\ 2 



(i = tw, dtw), 



(14) 
(15) 
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FIG. 5: Numerical solutions to the one-loop RG equations 
© for (aLJi/Ja = 0.2 and (b) J1/J2 = -0.2. We set A = 1 
[see Eq. (Ml)]. It is found that the three coupling constants 
Gi, Gtw, and Gdtw are most relevant and grow under the RG, 
with asymptotically a simple ratio Gi : Gt w : Gdtw =2:1: 
( — 1) or 2 : ( — 1) : 1 for J\ > and J\ < 0, respectively. In 
the plots, factors ±2 are muliplied to Gtw and Gdtw so that 
these ratios can be visually confirmed. 



Numerical solutions to the RG equations (16) are pre- 
sented in Fig. ||. For both signs of J\ , the three coupling 
constants Gi, Qfm and Gdtw finally grow to large values 
under the RG;t§E!l they asymptotically have the simple 
ratio Gi : G tw : G d tw = 2:1: (-1) or 2 : (-1) : 1 
for J 1 > and J\ < 0, respectively. Remarkably, Gi fi- 
nally grows with a positive sign for both signs of J\ . For 
J 1 < 0, in particular, it is initially negative but changes 
sign before starting to grow in the RG process. By con- 
trast, G tw retains the same sign as its initial value. The 
properties of the fixed points governed by large Gi(> 0), 
Gtw, and Gdtw are non-trivial. In fact, while the non- 
Abelian formalism allows us to derive the RG equations 
in a manifestly SU(2)-invariant form, it is often not very 
useful for discussing the physical roles of (marginally) rel- 
evant perturbations. In the next section, we proceed to 
the Abelian bosonization analysis to show that the pos- 
itive development of G\ induces a gapped state with a 
finite dimer order parameter -D123 7^ 0. 

As seen in Fig. ||(a) and (b), the coupling constants 
grow much more slowly for J± < than for J\ > 0. 
This implies that for J\ < 0, the energy gap associated 
with the dimer order should be much smaller and the 
spin correlation length £ should be much larger. In fact, 
as argued by Itoi and Qinj2il the correlation length be- 
comes of astronomical scale [e.g., £/a ~ e 83 ~ 10 36 for 
the case of Fig. Wb)]. Such a tiny gap or a large corre- 
lation length is very difficult to detect by any numerical 
investigation; the system effectively behaves like a gap- 
less system even when the system size is macroscopically 
large. We stress, however, that this insight is based on 
the perturbative RG analysis for small J1/J2 < 0, and 
it is possible that the energy gap grows to an observ- 
able magnitude as we increase \Ji\/J%. Our numerical 
result presented in Sec. Ill B indeed identifies a large but 



detectable correlation lengths around J1/J2 = —2. 



2. Abelian bosonization 



where A is a dimcnsionless constant of order unity. 
Using ihc-. oacrator product expansions in the WZW 
thggglgZrH the one- loop RG equations^ are derived 



Gbs — G 
Gi = Gi + Gt w — GtwGdtw, 



Gt w — ~-^Gb s Gt w - 
3 

Gdtw = ^^bsGdtw 



GiGtw 
-GiGt 



-G\G 



dtw 5 



(16a) 
(16b) 

(16c) 
(16d) 



where the dot indicates the derivative (Gi = dGi/dl) 
with respect to the change of the cutoff: a — > e dl a. See 
Appendix [A] for the derivation of Eq. ([l6]). We have 
ignored G2 since it does not affect the flow of the other 
coupling constants at the one-loop level. 



In this section, we use the Abelian bosonization for- 
malismcJ to discuss the physical roles of the marginally 
relevant perturbations Gi(> 0), Gt w , and Gdtw identi- 
fied in the non-Abelian analysis. Although the Abelian 
formalism obscures the SU(2) symmetry of the model, it 
has the advantage of simplifying identification of various 
orders with the pattern of locking of bosonic fields, as 
illustrated in Table |. 

Let us start from the two decoupled antiferromag- 
netic chains in the limit J1/J2 — ^i-Qr— 1 We summarize 
the Abelian bosonization descriptionEiEij of a single XXZ 
chain (0 < A < 1) , so t hat the same formulation can be 
used later in Sec. IV A . Each decoupled XXZ chain la- 
beled by n = 1, 2 is described by a Gaussian Hamiltonian 



Hn = I dx~ [K(d x e n ) 2 + K-\d x n f] (17) 
where the velocity v and the TLL parameter K are given 
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by 



2 arccos A 



1 — (l/V) arccos A 



(18) 



The bosonic fields <p n and 9 n satisfy the commutation 
relation 

[<p n {x),e nr {x')]=i5 nn ,Y{x-x l ), (19) 
where Y(x — x') is the step function 

(x < x'), 
Y(x-x') = { 1/2 \x = x'), (20) 

1 (a; > x'). 

The spin and (in-chain) dimer operators are expressed in 
terms of the bosonic fields as 



S. 



2j+n = —=d x 4>n{x n ) + (-l) 3 Ai cos[v27r^„(x n )] 
v 2n 



+Bi cos[v / 27T0„(x„)] + . . . I 



(-I) S 2 j+n ■ S 2j+n+2 = C sm(V2%(f> n ) + 



(21) 



(22) 
(23) 



where A u B , B l (Rcfs. |5lJ,[32j), and C (Ref. |53|) are non- 
universal constants which depend on A. 

We now focus on the case A = 1, at which K = 1. 
To treat the coupled chains, it is useful to introduce the 
bosonic fields for symmetric (+) and antisymmetric (— ) 
sectors: 

^± = -^±02), 0± = -j=(9i±02). (24) 

The three perturbations found to grow in the ..notu 
Abelian analysis have the following expressions Ho0 



B 

Oi=--\ cos(v / 470+) cos(V4^6L) 
2a z 

+ ^ [(^+) 2 - (d x 9 + ) 2 - (d x <t>-f + (d x 6-f] , 



(25a) 



+ 



2 . 



2n 



O, 



dtw 



^c 2 



(5 x + )sin(v47T(?!)_) 
+(d x (p-)sm(V^T(j> + ) 

{d x cj>+) sin(v / 47T(?!)-) 
-(d x (j>-) sin(v / 47r^+) . 



(25b) 



(25c) 



Furthermore, the 2 term, which is decoupled from the 
other terms in the RG equation (|l6|), has the expression 



02 



B 

- — \ cos(v / 47r0- ) cos(V47r6'_) 
2a z 

+ ^ [(d x <f>+) 2 + (d x e+) 2 - (a^_) 2 - (d x e.) 



2 ]- 

(26) 



The second lines of Eq. (25a) and Eq. (l2q) can be com- 
bined with the Gaussian Hamiltonians (uj\) of the decou- 
pled chains, leading to 

H = J dxJ^Y [ K -(9x9u) 2 + K-\d x tf> v f\ (27) 



u=± 



with 



#4 



V± = V 



it y + o(g 2 ,gI), 

l± 9l + {G\,Gl) 



(28) 



Using the new Gaussian Hamiltonian Hq, we can calcu- 
late the scaling dimension of the operators in Eq. (|25|). 

Specifically, the scaling dimension of e l ^*' i}± and e %y ^ e± 
is given by K ± and K ± 1 , respectively. In the non- Abelian 
analysis, we have seen that G\ grows to a positive value 
in the RG flow irrespective of the sign of J\. Assum- 
ing G\ > 0, we find that the pro duct of the two cosine 
operators in the first line of Eq. (25a) (with scaling di- 
mension 2 — G\) is the most relevant term among those 
in Eq. (Eq). This term locks the bosonic fields at 



(29) 



(V47r</) + , V4tt6L) = (0,0) or (tt.tt). 



These correspond respectively to finite positive or neg- 
ative value of the dimer order parameter D\ 2 s = (Si • 
S 2 ) — (S 2 ■ S3}, since the (inter-chain) dimer operator is 
expressed as 

S 2 j+i • S 2 j +2 — S 2 j+2 ■ S 2 j+3 = 2a 2 Ni ■ N 2 + . . . 

w 2B 2 cos(v / 47r6'_) + A\ cos(v / 47r0+) + cos(%/47r</>_) 

(30) 

In the last expression, the first term and the rest come 
from the xy and z components of the spins, respectively. 
For the locking in Eq. (|2^), these components acquire 
both positive or both negative expectation values, in 
agreement with Fig. |^(c) and with the SU(2) symme- 
try of the model. It is worth noting that the locking 
positions of the two degenerate ground states in Eq. ( |29| ) 
are independent of the sign of J\ in the isotropic case 
A = 1. The second most relevant terms in Eq. ( j25| ) are 
(d x Q+) sin(A/4~7r#-) and (d x (j>-) sin(V^4~7r0+ ) wit h scaling 
dimension 2 — G\j2. As explained in Sec. IV A, the for- 
mer has the effect, of inducing the incommensurability in 
spin correlations E Since a finite energy gap opens due to 
Gi > in the dimer phases, the incommensurate spin 
correlations are expected to remain short-ranged. 
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B. Numerical results and physical properties of 
dimer phases 

In this section, we present numerical results on the 
model (|l|) in the isotropic case A = 1, and discuss physi- 
cal properties of the dimer phases for different signs of J\ . 
In agreement with the field-theoretical results reviewed 
in the previous section, we find that the dimer order pa- 
rameter -D123 becomes finite for both signs of J\ , and that 
there are doubly degenerate ground states with positive 
and negative Di 23 . While we propose different physical 
pictures fo r the d imer orders in the J\ > and Ji < 
cases (Sec. IIIB 1 ), we also discus s a hidden order com- 
mon to the two cases (Sec. IIIB 3). In the following, our 
numerical results (based on iTEBD with \ = 300) are 
presented for the ground state with -D123 < 0. 



1. Local spin correlations 

In Fig. ||(a), we plot nearest-neighbor spin correlations 
(Se ■ Se+i) (with £ = 1, 2) and the dimer order parameter 
D 123 = (Si ■ S 2 ) - (S 2 ■ S 3 ) for -3 < J1/J2 < 3. While 
-D123 7^ can be confirmed for both J\ > and J\ < 0, 
a notable difference between the two cases can be found 
in the signs of local spin correlations. 

For J 1 > 0, one of the following inequalities is always 
satisfied: 

(Si -S 2 ) < -(S 2 -S 3 ) <0 (0<Ji/J 2 <2), (31a) 
(Si ■ S 2 ) < (S 2 ■ S 3 ) < (2 < J1/J2 < 4.15). (31b) 

Namely, the system has a strong antiferromagnetic cor- 
relation on the bond (1,2) and a weaker correlation on 
(2,3). In this case, it is natural to assume that singlet 
dimers are formed on the bonds (2j + 1, 2j + 2) (j € Z), 
and are weakly correlated with each other, as schemat- 
ically shown in Fig. 0(a). Hence we call this phase the 
singlet dimer phase. In particular, the ground state is 
exactly given by a direcLproduct of singlet dimers at the 
Majumdar-Ghosh pointH Ji/J 2 = 2. In Fig. |(a) we find 
that the weaker correlation (S 2 ■ S3) changes the sign at 
this point. 

By contrast, the following inequality is found to be 
satisfied when —4 < J\j J 2 < 0: 



Haldane dimer 



< {Si ■ S 2 ) < (S 2 ■ S 



■V- 



(32) 



Namely, strong and weak ferromagnetic nearest-neighbor 
correlations alternate along the chain. This observation 
led us to propose that there should be emergent spin-1 
degrees of freedom on the bonds (2j + 2, 2j + 3) (j € Z) 
that have stronger ferromagnetic correlation, as depicted 
by ellipses in Fig. 0(b). Since the total wave function is 
a spin singlet, such spin-l's are expected to form a va=. 
lence bond solid stateo as in the spin-1 Haldane chain.Ej 
Namely, from each encircled bond in Fig. ^|(b), two va- 
lence bonds emanate, one to the left and one to the right; 
the total wave function is obtained by superposing such 
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FIG. 6: (Color online) (a) Nearest-neighbor spin correlations 
(Sj ■ Sj+i) and the dimer order parameter D123, (b) the spin 
correlation length £, and (c) the string correlation (Q) with 
t — 2 and r = 50, as a function of Ji / J2 in the isotropic case 
A = 1. In panel (a), D123 multiplied by 40 is also plotted 
for J1/J2 < (filled circular symbols). In panel (b), £ is too 
small to determine around J1/J2 = 2. 



valence bond covering states. We thus call the dimer 
phase with J\ < the Haldane dimer phase. The emer- 
gence of the Haldane chain physics in this phase is also 
supported by th e prese nce of a hidden non-local order 
analyzed in Sec. Ill B 3| . 



In Sec. Ill A, it was argued that the marginal perturba- 
tion G\ , which induces the dimer order, grows very slowly 
under the RG for J\ < and that the energy gap associ- 
ated with the dimer order can be extremely small.cil The 
result of Fig. ^(a) indicates that the dimer order param- 
eter D123 grows to a numerically detectable magnitude 
for intermediate values of \Ji\/J 2 (~ 2), although the ob- 
tained values are much much smaller compared to the 
Ji > case (by a factor of around 1/40). The weakness 
of the effect of J\ in inducing the dimer order and the 
associated energy gap for J\ < is also seen in the spin 
correlation length discussed next. 
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(a) Singlet dimer state (J^O) 





2 ; [4 














5 j 
) 1 



(b) Haldane dimer state (J^O) 




FIG. 7: (Color online) Sketches of (a) the singlet dimer state 
and (b) the Haldane dimer state. The thick lines indicate 
valence bonds. In (b), the encircled bonds indicate emergent 
spin-l's. From each of them, a valence bond emanate to each 
left and right; the wave function is given by a superposition 
of such valence bond covering states. Vertical cuts (dashed 
lines) are introduced to probe a hidden order; the number of 
valence bonds crossing with each cut is shown in the square. 
The alternation of odd and even numbers is found in both the 
states. 



2. Spin correlation length 

We determine the spin correlation length £ in the dimer 
phases by using the method of Ref. ||. Except at the 
Lifshitz point J1/J2 = 2, the spin correlation function is 
expected to behave at long distances aa§E3 



{Si ■ Si+r) w Acos(Qr)r-h- r/(i . 



(33) 



In the incommensurate regions —4 < J1/J2 < and < 
J1/J2 < 2, the pitch angle Q chang es cont inuously from 
to 7T, as will be discussed in Sec. IV B 5| (sec Fig. 16). 
For 2 < J1/J2 < 4.15, Q is fixed at Q = it. To determine 
£, we plot r 1 ' 2 e r '^{S\ ■ S\ +r ) as a function of r, and 
tune £ such that the amplitude of oscillations becomes 
as constant as possible, as illustrated in Fig. [|. While 
the coefficient A in Eq. ([}3|) is given by the oscillation 
amplitude in Fig. ^, it is not simple to determine Q which 
can fit these very rapid oscillations; instead it will be 
determined by calculating the spin structure factor in 

Fig. m 

The calculated £ is plotted in Fig. ||(b). Tbfi data for 
Ji > are broadly in agreement with Ref. ||.lZI We find 
that the values of £ are much larger for J\ < than for 
J 1 > 0, as anticipated from the magnitudes of the dimer 
order parameter in Fig. [|(a). 
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FIG. 8: Determination of the spin correlation length £, illus- 
trated for J1/J2 — —1.8. Assuming the asymptotic behavior 
@, we plot the function r 1/2 e r/4 (5i ■ Si +r ), and tune ^ 
such that the oscillation width of th*« function becomes as 
constant as possible as a function of ru While the oscillations 
arise from the cosine factor in (^), it is difficult to extract 
the pitch a ngle Q from this figure; instead, calculations in 
Sec. 



IVBj give Q/(2tt) w 0.235. 



We use the above numerical data of the spin correla- 
tion length £ to infer the magnitude of the spin gap A s 
for Ji < 0. In general the spin gap A s should be inversely 
proportional to £, with the proportionality constant be- 
ing the spin velocity. From the data of Ref. || for J± > 0, 
we extract an approximate relation (A s /J2)£ ~ 2. Ap- 
plying the same relation to the J\ < case, we estimate 
the spin gap A s around J1/J2 = —2 to be roughly equal 
to O.O6J2. We note that this should be considered as a 
crude order of magnitude estimate. 



3. Hidden order 

The singlet and Haldane dimer phases have different 
(local) features of short-range correlations as expressed 
in Eqs. (|3l]) and (^2|). In spite of this local difference, 
the two phases in fact share a common non-local order, 
as we now explain. Let us count the number of valence 
bonds crossing the vertical cuts (dashed lines) depicted in 
Fig. @. We find that even and odd numbers alternate in 
the same way in the two phases, when we take the ground 
state with D123 < 0. The existence of such a hidden non- 
local order can be probed pupiprically by calculating the 
string correlation function&tS 



o z st A 



J. + 2r) 



(S! + Si +1 )exp 



£+2r-l 
m=e+2 



x (5, 



i+2r 



S: 



£+2r- 



(34) 
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Og tr (2,2+2r): J 1 /J 2 = 
J 1 /J 2 =- 
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(a) J 1 /J 2 = 1.4 



30 r 40 



FIG. 9: (Color online) String correlation function ( p4| ) for 
J1/J2 = 1.4 and —2.0 in the isotropic case A = 1. For 
both values of J1/J2, O z (2, 2 + 2r) remains finite in the long- 
distance limit while O z (l, 1 + 2r) decays to zero. 



The intuition behind this expression is as follows. Con- 



sider a pair of spins S t 



e+2j 



on the bond 



2j, I + 2j + 1) (j € Z), which the string correlation func- 
tion (|34|) consists of. If an odd number of valence bonds 
cross any cut placed between the neighboring pairs, then 
the pattern of Sf +2 j + 5 , | +2j+1 = —1, 0, +1 shows a hid- 
den antifcrromagnetic order, namely, alternation of +1 
and —1 after removing all 0's (see figures in Refs. ^l| and 
|63| ) . The correlation function ( |34| ) detects this hidden or- 
der and takes a non-vanishing value in the long-distance 
limit r — > 00. 

Figure ^ presents the numerical data of the string cor- 
relation functions ( |3~4| ) calculated with different starting 
points £ = 1,2 for the ground state with -D123 < 0. We 
find that for both signs of Ji, O z (2, 2 + 2r) remains finite 
in the long-distance limit while O z (l,l + 2r) decays to 
zero, in agreement with the even-odd structure in Fig. |?]. 
We note that this behavior is also consistent -ssnth the 
bosonized expressions of the string correlationsEj 



O 



str (l,l + 2r) 



O s z tl -(2,2- 



■2r) 



(cos[v / ^0+(x)]cos[v /: ^> + (y)]), (35) 
(sin[ v / i<M20]sin[ v / ^<My)]), (36) 



(with x and y being the two endpoints of the string) and 
the field locking position \Z4n(p + = it for the ground state 
with L»i2 3 < [see Eq. {^}. The J\j J 2 -dependence 
of O z (2,2 + 2r) for a long distance r = 50 is shown in 
Fig. [j](c). Although the dimer order parameter shows a 
large difference in magnitude between the J\ > and 
Ji < cases, the values of the string correlation are 
rather comparable between the two cases. 

Another way of probing the hidden order i&-,to find 
the degeneracy in the entanglement spectrumEj Using 
the Schmidt coefficients {A mQ } calculated in iTEBD, we 
plot the spectra {— 21nA mQ } in Fig. [l(]. Here the spectra 
arc classified by the z-componcnt magnetization to in the 
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FIG. 10: (Color online) Entanglement spectra {—2 In X ma } for 
(a) J1/J2 = 1.4 and (b) J1/J2 = —2. The left and right panels 
are for the bipartition of the system at the bonds (1,2) and 
(2, 3), respectively, m refers to the magnetization in the right 
half of the system. The lower entanglement level corresponds 
to the more important weight in the total state. We note that 
as an example, the exact singlet dimer ground state of Ref. [| 
shows — 21nA ± i 1 — In 2 and — 21nAo,i = for the two types 
of bipartition (with all the other levels at infinity) . 



right half of the system (this classification is done in the 
process of our calculations to exploit the U(l) spin rota- 
tional symmetry for better efficiency) . For the bipartition 
of the system at the bond (1, 2) (left panels), we find that 
the entanglement levels appear only for half-integer m, 
and are all doubly degenerate due to the left-right sym- 
metry around m = 0. By contrast, for the bipartition at 
(2, 3) (right panels), the entanglement levels appear only 
for integer m, and non-degenerate levels are found for 
to = O.Ej These features are found commonly for both 
signs of Ji, and are consistent with the even-odd struc- 
ture in Fig. ^. 

In Fig. 0, we depicted short-range valence bonds only. 
However, the even-odd structure we discussed can be also 
defined in the presence of longer-range valence bonds. As 
the correlation length becomes longer, the weights of such 



12 



D m > 0, 

= (°.°) 

mm 



AAA 



j[u 2 





FIG. 11: Expected phase diagram of the zigzag ladder 
model with alternating nearest-neighbor couplings Ji and Ji . 
IJ1I/J2 and IJ1I/J2 are assumed to be small. The solid di- 
agonal line Ji = J[ corresponds to the original model (|l]) 
(with A = 1), and represents the first-order phase transition 
line in the current model. The vertical and horizontal dashed 
lines correspond to a usual ladder model (no phase transition 
on these lines). In four insets of zigzag ladders, thick lines 
indicate valence bonds, and ovals indicate the formation of 
effective spin-l's. 



longer-range valence bonds in the wave function would 
gradually grow while retaining the even-odd structure.EH 
We expect that through this process, the Haldane dimcr 
state of Fig. 0(b) smoothly changes into the exact res- 
onating valence bond ground state at J1/J2 = — 4, in 
which vaLeace bonds are uniformly distributed over all 
distances £3 



model (the lower half of the vertical axis of Fig. |lT|) by 
gradually switching off the J\ coupling. It is also possible 
to adiabatically change the ground state from the Hal- 
dane dimer state to the singlet dimer state (both with 
D123 < 0) by moving counterclockwise around the origin 
in Fig. [n], although the wave function may considerably 
change in this process. In the zigzag ladder model with 
Ji = J{ (diagonal line), however, the singlet and Haldane 
dimer phases are separated by the origin (open circle in 
Fig. |ll]), at which the two chains are decoupled. We note 
that only on the J\ = J[ line in Fig. 11, the model has 
the symmetry with respect to the translation Si — > Se+i, 
and the dimer order appears by spontaneously breaking 
this symmetry. It would thus be interesting to investi- 
gate under what kind of translationally symmetric per- 
turbation the Haldane and singlet dimer phases can be 
adiabatically connected to each other while retaining the 
double degeneracy below a finite excitation gap. 



IV. EASY-PLANE CASE < A < 1 

In this section, wc consider the mod el fljl ) in the easy- 
plane case < A < 1. In Sec. [IV A| , we present 
the Abelian bosonization formulation of the model for 
I J\ I / J2 *C 1 and explain how various phases in Fig. [t] are 
described in this framework. In particular, we review the 
effective theory for the gapless chiral phasetl and, follow- 
ing Ref. [54], discuss its instability towards gapped chiral 
phas es du e to a symmetry-allowed perturbation. Scc- 
1VB presents our numerical results. We compute 



tion 

the spin correlation functions in the gapless chiral phase 
and determine the phase boundaries to the gapped chiral 
phases. 



4- Adiabatic connectivity to a ladder model 



In order to gain further intuition about the two dimer 
phases, it is useful to introduce explicit bond alternation 
of the J 1 couplings in the Hamiltonian (Q) (with A = 1). 
Namely, we place inequivalcnt couplings J\ and J[ on 
the bonds (2j + 1, 2j + 2) and (2j + 2, 2j + 3) (j g Z), 
respectively Figure 11 displays an expected phase di- 
agram for small > hi J2 and J[/ J2- This phase diagram 
can be obtainedrir 3 ! by noticing that in the non- Abelian 
bosonization framework, the bond alternation induces 
the relevant term (Ji — J[)Ni ■ iV 2 with scaling dimen- 
sion 1 in the Hamiltonian, which leads to the ground 
state where D123 ~ (iVi • iV 2 ) acquires a finite average 
with the same sign as that of J{ — J±. The limit Ji — > 
or J[ —> (the vertical or horizontal axis of Fig. 0) 
corresponds to a spin ladder model, for which it is estab- 
lished that the rung singlet and Haldane phases appear 
for antifcrro maMiCj iic and ferromagnetic rung couplings, 
rcspectively.cZl£ili£3'EJ Therefore, we expect that the Hal- 
dane dimer state with -D123 < in Fig. 0(b) should be 
adiabatically connected to the Haldane state of a ladder 



A. Bosonization analyses 

We consider the easy-plane XXZ Hamiltonian ([l]) in 
the reg ime \ J\ \/J% *C 1. Using the formulation described 
in Sec. Ill A 2, we obtain the effective Hamiltonian 



i/=± 



— 71 cos(V47T(/> + ) cos(V47r6L) ^37^ 
+ 7tw(dz^)sin(v / 4^6L) 

+ 7^(^_)Bm(V4^+) + •••}. 

The first line represents the Gaussian Hamiltonian while 
the other lines represent pertuijb|&lioi|is which can become 
relevant in the easy-plane case.Bn3a As seen in Eq. (|2^), 
the 71 term is related to the G\ term in the non- Abelian 
bosonization, while 7 tw and 7^ correspond to Gt w - The 
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coupling constants are obtained in lowest order in J\ as 



K[ 1 T 7T~ I ,«± 



2irv 



Ti 



7tw 



7tv 



l± ~tr^ (38) 



^-JiA^, (39) 



where K and u are given by Eq. (^8|). We have discussed 
in Sec. Ill that, in the isotropic case A = 1, 71 grows to 
large positive values for both signs of Ji under the RG, 
and induces the singlet and Haldane dimer phases for 
J\ > and Ji < 0, respectively. Below we explain how 
other phases in Fig. [l] are described using the effective 
Hamiltonian ( |37| ) . The results are summarized in Table |. 



1. Even-parity Aimer phase 

If J 1 < 0, the coupling constant 71 is negative at the 
bare level. Suppose that this term grows, keeping the 
negative sign under the RG. Then the bosonic fields are 
locked at 



(V47T0+, V47r6>_) = (0,tt) or (vr,0). 



(40) 



In either case, it follows from Eq. (|30|) that the xy and 
z components of the dimer order parameter, D^% 3 and 
DI23, become finite and have mutually opposite signs 
(D123DI23 < 0). This situation corresponds to the 
even-parity dimer phase appearing at strong easy-plane 
anisotropy (A < 0.6); see Fig. ||(b). 



2. Gapless chiral phase 

As shown by Nersesyan et aL,i the gapless chiral phase 
appears when 7 tw grows under the RG. To discuss the 
effect of the 7^ term, it is useful to perform the mean- 
field decoupling 



(41) 



(d x 6 + )sm(V4:Tr6-) + (d x 9 + )(sm(V<liT9-)) 



Then the Hamiltonian (|37| ) separates into "+" and "— " 
sectors: 



H = H< 



with 



dx— 



K+{d x e + f + K+ 1 (d x( f )+ ) 2 

J dx{^- [K_(d x 9_) 2 + KZ\d x ^)- 
+ 7tw(d*0+}sin(V4^0_)}. 
Here we have introduced 

7t w (sin(\/47r0-)) 



qx, 



(42) 
(43) 

(44) 
(45) 



While H + is a Gaussian Hamiltonian of free bosons 
(<t>+, 0+), H_ is a sine-Gordon Hamiltonian in which the 
relevant sine potential generates a finite energy gap for 
the 6L field. Since (dj + ) = from H+, (d x 0+) = q. 
The coefficient of the sine potential in H- is thus given 
by 7tw<7, and the field 0_ is locked at distinct positions 
depending on the sign of this coefficient: 



(46) 



(V4tt6L) = --sgn(7 tw g). 



Correspondingly, the sine term acquires a finite expecta- 
tion value: 



(47) 



(sin(V47r0_)) = -ci sgn(7 tw g), 



where c\ is a positive constant__Equations ([45]) and ( |47| ) 
can be solved self-consistentlyLJ by inserting the exact 
solution of the sine-Gordon model into Eq. (47), yielding 
two solutions, one positive and one negative q. It should 
be understood that the mean-field parameters c\ and q 
used in the following calculation of correlation functions 
are determined selfconsistcnly. 

First, the non-vanishing value of the mean-field pa- 
rameter in Eq. (B7|) directly leads to a finite vector chiral 
order parameter (0): 



K i,i+i = -Bo( sin (v47r6L)) = B 2 Cl sgn( 7tw a). (48) 

Therefore the two mean-field solutions correspond to the 
ground states with positive and negative k| Let us 
take the ground state with n\ e+1 > (i.e., 7t w g > 0) and 
discuss the expressions of the spin operators. We focus 
on gapless degrees of freedom, and ignore the fluctuations 
of 9- around its average (Eq). Then we find 



27r(9i i2 = Vn{0+ ±9-) = V¥0+ + y/nqxi^ T 7, (49) 



which arc combined into 



2ir9 n (x n ) = \fKd+{x n ) + \pKqx n + — 



3 

n ~2 



(50) 



The in-plane component of the spins are then expressed 
as 



S, 



2j+n 



B (-iy e iV ^ 8n ^ 



Bo cxp 



V^(0+ + qx n ) + - [2j + 



(51) 



Introducing I = 2j+n and x(t) = x n (j) = (a/2) (£-3/2), 
we obtain 



S /, 



B e l[ ^ §+{x}+Q< - £ ~ 3/2) K 



with 



Q 



V^qa 



(52) 



(53) 
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TABLE I: Summary of the Abelian bosonization description of the phases for small \Ji\/ J 2 - Both the easy- plane (Sec. |Fv|) and 
easy-axis (Sec. |v|) cases are presented. We note that the (chiral) even-parity dimer and chiral Neel phases appear for rather 
large \Ji\/J 2 in Fig. | although their essential features can be captured in the Abelian bosonization framework. 

Phase Relevant perturbations Field-locking positions Order parameters 

Singlet /Haldane dimer 71 > (y/4n<j>+,y/4n9-) = (0,0), (7r,7r) D*% 3 Df 23 > 

Even-parity dimer 71 < (V4^(t>+,V4^0-) = (0,7r), (tt,0) Df% 3 Df 23 < 

Gapless chiral 7 t w ~ J\ 7^ = — ^sgn(Ji(d x 9+}) k\ 2 7^ 

Chiral singlet/Haldane dimer 7tw ~ Ji, 7nd < 0, 71 > (V^<t>+, \/47r0_) = \ ^ 0,± ^ ^ ( - 0,0 ^ k\ 2 0, Df% 3 Df 23 > 

I (7T,±f) -> (7T,±7T) 

Chiral even-parity dimer 7 tw < 0, 7 nd < 0, 71 < (V4jv4>+, V^kQ-) = \ ^ 2 «*2 / 0, D*% 3 Df 2S < 

[(7T,±f) (?r,0) 

Chiral Neel 7tw < 0, 7nd > 0, 7^ + V^</>+ = ±tt/2, = ±tt/2 k\ 2 / 0, (Sf ) oc (-1)* 



Uudd 7bs < V^-l = 0, 7T, V2^02 = 0, 7T (S 2 Z 3 + 1 ) = ±{SQ OC (-1) J 

Partially polarized j' tv < v 7 ^- = f sgn((9 ;c + }) (S|) / 



As for the z component of the spins, we simply ignore 
the </>_ part of the expression: 



Sf 



47T 



d x (, 



Spin correlation functions are then calculated a; 

e -iQ{l'-t) 
{S+S-,) = A 2K+) + . . 



(SIS!,) 



K+ 
2tt 2 \£'-£\ : 



(54) 

I 

(55) 
(56) 



with A = B 2 2 1/( - 2K +' ) . The finite vector chiral order pa- 
rameter Kg i+1 in Eq. (^) and the quasi-long-range in- 
plane spiral correlation with an incommensurate pitch 
angle Q in Eq. ( pq ) are two major features of the gapless 
chiral phase. 



3. Gapped chiral phases 

Following Lecheminant et a/.,@ we consider the follow- 
ing symmetry-allowed perturbation to the effective the- 
ory of the gapless chiral phase: 



7nd 



J dx cos(2' V / 47r0+), 



(57) 



with which the "+" sector of the Hamiltonian becomes 
a sine-Gordon model. The scaling dimension of this per- 
turbation is 4A' + . If the 7„d term becomes relevant 
(4JC+ < 2), a Berezinskii-Kosterlitz-Thouless (BKT) 
transition takes place and as a result, the bosonic field </>+ 
is locked at distinct positions dependent on the sign of 
7nd- This leads to gapped chiral phases in which the chi- 
ral order coexist with either the dimer or the Neel order, 
depending on the sign of j n( \. 



First, when 7 nc i < 0, \ / 4tt(/) + is locked at 



which produces a finite value of the ^-component of the 
dimer order parameter, Df 2 3, as seen in Eq. (|30|). We 
have thus obtained the "chiral dimer phase," in which 
the vector chiral and dimer orders coexist.E3 Once 0+ is 
locked as in Eq. dSq) , the locking position of y/4ir9^ is 
affected by the 71 term in Eq. pTl) and changed from 
±7r/2 of the gapless chiral phase [Eq. ((itf)], so that the 
xy-component of the dimer order parameter, -D^, also 
becomes finite, in agreement with Fig. ^. Specifically, for 
positive 71 , the field- locking positions of the four degen- 
erate ground states change smoothly with the strength 
of 71 as 



(0,±f) 



(0,0), 

(-7T, ±7t) 



(7T,7r) 



(59) 

finally resulting in the two degenerate ground states of 
either the Haldane dimer or the singlet dimer phase as 
specified by Eq. (p9|). For negative 71, the field- locking 
positions change as 



M), 



(V4tt0 + , V4tt6L) 




(60) 

resulting in the two degenerate ground states of the even- 
parity dimer phase as indicated by Eq. (|40|). 
Second, when 7 n d > 0, 4>+ is locked at 



(61) 



This yields a finite Neel order paramter along the z di- 
rection, as we explain below. From Eqs. ( ^l[ ) and (|24|), 
the Sg operator has the staggered component 



47T 



(62) 



47r</>+ = or 7T, 



which, at first sight, looks insensitive to the locking of 
(j> + . However, after the locking (|6l|), the 7^ term in 
(58) Eq. ([37]) reduces to the operator ±j" w d x cf>-, which can 
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FIG. 12: (Color online) In-plane spin correlation function 
\{Sf S^, r }\ for fixed A = 0.8 and various values of J1/J2 in 
the gapless chiral phase. Logarithmic scales are used in both 
axes. 



be absorbed into the Gaussian part of in Eq. (|44|) 
by redefining <f>_ (so that d x 4>_ is shifted by a constant). 
Consequently, Eq. (|6^) acquires a nonvanishing expecta- 
tion value 



(-iY(sf) = 



K+a 



7t ';(sin(V47r</» + )) 



(63) 



We have therefore obtained the "chiral Neel phase," in 
which the vector chiral and Neel orders coexist. 

At the BKT transition point K + = 1/2, the sine- 
Gordon theory for the "+" sector predicts the appearance 
of a multiplic aiixe .logarithmic correction to the correla- 
tion functionsMyOO 




-0.25 0.25 

q / (2k) 



FIG. 13: (Color online) Equal-time spin structure factors 
[Eq. @], (a) S + ~(q) and (b) S zz (q), in the gapless chiral 
phase. Calculations were done for the same parameter points 
as in Fig. and we set L = 100. 



<# S *> = A & \l>-l\ lnV2(|£ ' " el/a) + ' ' ' ' (64) 

This logarithmic correction is utilized to locate the BKT 
phase transition point numerically in the next section. 



B. Numerical results 

In this section, we present our numerical iTEBD re- 
sults (with the Schmidt rank x = 300) on the spin corre- 
lation functions in the easy-plane case < A < 1. 



1. Spin correlations in the gapless chiral phase 

We first discuss the numerical results for the gapless 
chiral phase, where we choose the ground state with 
K12 > 0. Figure [l2] shows the in-plane spin correla- 
tion function | {SfS^ +r ) | at A = 0.8 for various val- 
ues of J1/J2 in the gapless chiral phase. The data for 
\Ji \/ Ji <L 2 follow straight lines in logarithmic scales, in 



agreement with the power-law behavior in Eq. (|55|). By 
contrast, the data for J1/J2 = —1.6 and —1.2 show some 
oscillations at short distances although the overall behav- 
iors are linear as expected from Eq. ([55]) (we suspect that 
the downward bending at large r for J\j J2 = — 1.2 is due 
to a finite Schmidt rank \ = 300, and is not a genuine 
behavior). 

The origin of the oscillations can be found in the spin 
structure factors shown in Fig. |l3|. For L consecutive 
spins at the sites I = 1, 2, ... ,L in a translationally in- 
variant infinite system treated by iTEBD, we introduce 



(65) 



and define the equal-time spin structure factors as 

S^(q) = (S°S f l q ) with (a,P) = (+,-), (z,z). (66) 

In Fig. fj"3|(a), S~* (q) shows sharp peaks at incommensu- 
rate wave number q = Q > 0, which become sharper and 
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FIG. 15: (Color online) In-plane spin correlation, calculated 
for fixed A = 0.8 and various values of J1/J2 around the 
transition points shown in Fig. ^| The symbols "C" and "N" 
indicate our estimates of the transition points (with a pre- 
cision of 0.005) for the onsets of the vector chiral and Neel 
orders, respectively. 



it shows finite peaks at incommensurate q. Although 
the explanation of these peaks is beyond the scope of 
the effective theory, their occurrence is rather natural for 
A = 0.8, since the xy and z components should show 
similar behaviors as the system approaches the isotropic 
limit A = 1. 



FIG. 14: (Color online) In-plane spin correlation function, cal- 
culated for fixed J1/J2 = — 2 and various values of A around 
the transition points shown in Fig. [2]. The symbols "C" and 
"D" indicate our estimates of the transition points (with a 
precision of 0.005) for the onsets of the vector chiral and dimer 
orders, respectively. Logarithmic scale is used for the horizon- 
tal axis. At the BKT transition related to the onset of the 
dimer order, the plotted function is expected to become linear 
in the long-distance limit, which we use to determine the "D" 
points. In panel (b), all the curves are slightly bent downward 
around r = 70 (broken vertical line) due to the finiteness of 
the Schmidt rank x(= 300) in iTEBD, so we use the range 
r < 70 for our analysis. 



higher for large \Ji\/J-z- This feature is consistent with 
Eq. ([35]) , provided that K + becomes larger with increas- 
ing IJ1I/J2; see Eq. (|38|). These peaks are expected to 
diverge as L — > 00 in the gapless chiral phase. For small 
IJ1I/J2, a second peak around q = — Q < develops, 
which indicates the ellipticity of the spiral correlations 
and is the origin of the oscillating behavior in Fig. [l|. 
The appearance of the second peak can be understood 
by observing that S A (q) should gradually become left- 
right sympetric as the vector chiral order parameter k z 
decreases.c3 In Fig. |l3|(b), S zz (q) shows linear behaviors 
around q = as expected from the Fourier transform of 
Eq. @: S zz {q) = K+\q\/2ir for |g| < 1. In addition, 



2. Transitions to the gapped chiral phases 

Next we analyze how the spin correlation changes at 
the transition from the gapless chiral phase to the gapped 
dimer or Neel phase. The existence of the intermediate 
gapped chiral phases where two kinds of orders coexist 
is anticipated from the analyses of the order parameters 
and entanglement entropy in Figs. p| and || and from the 
bosonization analysis of Sec. IV A 3|. The in-plane spin 



correlation function is expected to show a multiplicative 
logarithmic correction in Eq. ( |6~4| ) at the BKT transition 
point from the gapless to gapped chiral phases. There- 
fore, in Figs. and P, we plot r 2 |(5+ S{ 



+r/ 



which is 

expected to become a linear function of In r at the BKT 
transition point. In Figs. 0(a) and (b), the symbols 
"C" indicate the Ising transition points (determined in 
Fig. ^ at which the inversion symmetry is spontaneously 
broken and the vector chiral order appears. Finding the 
linear behavior of the plotted functions, we determine the 
BKT transition points as indicated by the symbols "D" . 
Narrow but finite ranges of intermediate phases between 
"C" and "D" are found in the intervals 0.61 < A < 0.63 
and 0.91 < A < 0.92, which we identify with the "chi- 
ral (even-parity and Haldanc) dimer phases." Similarly, 
we determine the range of the "chiral Neel phase" in 
Fig. Oil In this way, we have determined the "x" sym- 
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FIG. 16: (Color online) Pitch angle Q as a function of J1/J2 
for different values of A. This angle is determined by finding 
the peak in the in-plane structure factor S + ~(q) as shown in 
Fig. [l3](a). The classical value Q = arccos(— jjj), which is 
independent of A, is plotted together for comparison. 



V. EASY- AXIS CASE A > 1 

To complete our analysis of the XXZ chain model (jl|) 
with J\ < and J 2 > 0, let us shortly discuss the caae 
with easy-axis anisotmpy A > 1. In this case, IgarashO 
and Tonegawa et al.i3 have found the following three 
phases. For J\ j J-x < —4, the ground state is fully po- 
larized (ferromagnetic) along the z direction. For small 
\J\\I J-2 and large A, the ground state is antiferromag- 
netic, having a period-4 structure WW . . . (uudd). Be- 
tween the fully polarized and uudd phases intervenes the 
partially polarized phase, in which the spontaneous fer- 
romagnetic moment along the z direction changes contin- 
uously as a function of J1/J2 and A. We note that the 
found in the model with antifcrro- 
Here we describe the uudd and par- 
tially polarized phases in terms of the Abelian bosoniza- 
tion formulation for | Ji|/J 2 <C 1 and < A - 1 « 1. 



A. uudd phase 



uudd phase was a 
magnetic J\o, > 



bols in Fig. 0. Since the method of determining the BKT 
point from the logarithmic correction to spin correlation, 
as employed here, has not been discussed in literature 
(as far as we know), we demonstrate its validity using a 
simpler example in Appendix M. 



Pitch 



Finally, we determine the pitch angle Q of the in- 
commensurate spin correlations in the vector chiral and 
gapped phases. It is determined from the maximum 
position of the in-plane structure factor (q) (as in 
Ref. |74|). The data of Q so obtained as a function of 
J1/J2 are shown for different values of A in Fig. [T|. 
The Lifshitz points, at which the in-plane spin correla- 
tion function changes its character from incommensurate 
to commensurate (Q = or it), occur inside the singlet 
dimer phase for Ji > and inside the even-parity dimer 
or Neel phase for J\ < 0. For Ji > and all values of A, 
the determined Lifshitz points are very close to the point 
J1/J2 = 2 with the exact singlet dimer ground states. 
According to the argument of Ref. the Lifshitz points 
should be in fact located exactly at J\/ J2 = 2. The small 
discrepancy comes from the difference in the definition 
of Q; in Ref. it is defined in terms of the asymptotic 
behavior of the correlation function in the long-distance 
limit. For J\ < 0, the determined Lifshitz line is drawn 
by broken Uap, in Fig. [|; it starts from the highly degen- 
erate pointHo (Ji/J 2 , A) = (—4, 1) and ends near the 
point (Ji/J2,A) = (—2,0) with the exact even-parity 
dimer ground states. 



We start from the decoupled isotropic Heisenberg 
chains with J2 > 0. The in-chain easy-axis anisotropy 
J 2 (A - 1) J2 jin Sl j+n Sl j+n+? (with A > 1) adds to the 
Hamiltonian the backscattering terms 



7bs 



cos(2V27T0i) + cos(2V27rci2) 



(67) 



with 7bs < 0. If this term grows dominantly under the 
RG, the fields are locked at 



(V27T0i,V27r02) = (0,0), (0,tt), (tt,0), or(7r,7r). (68) 

These four-fold degenerate ground states correspond to 
the period-4 uudd structures with 

(si j+l ) = c 2 (-iy, (5| i+2 ) = ±c2(-iy, (69) 

where C2 is a non-zero constant [see Eq. (plj)]. 

B. Partially polarized phase 

The partially polarized phase found numcricallyEl can 
be understood from the m-cap-field treatment of the op- 
erator ( d x <t>+ ) sin(y / 47r^)_)J^3 which is contained in Ot w 
in Eq,— I j25h|) . Here we review the formulation of Zarca 
et al.f.^ and then discuss the behaviors of correlation 
functions., which were not discussed in detail in previous 
studiesMBa 

We start from the effective Hamiltoniana 



+ 7 t ' w (d x< £ + )sin(\/4^_)} 



(70) 
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with 7j w < 0. T he mea n-field decoupling similar to the 
one used in Sec. IV A 2 yields the effective Hamiltonian 
H = H-L. + H_ , where 



H, = I dx^± 



K+(d x e + ) 2 +K^(d x ct> + ) 2 \, (71) 
H_= f dx{^ [K-{d x 0-f + KZ 1 {d,A-) 2 } 
+ 7tw/ iSm (v / 47r0-)}- 
Here we have introduced 



(72) 



<f>+(x) = </>+(x) - fix, 



with 



fi = = {d x <t>+). (74) 



There are two self-consistent solutions: ft = \fi\. 
A non- vanishing fi directly leads to the spontaneous mag- 
netization 



(S|)= " <^ + >= " M = M. (75) 



Furthermore, the sine potential in H_ locks the bosonic 
field at 



(76) 



To see the physical consequence of the field locking 
in Eq. (76), we discuss spin correlation functions in the 
ground state. The transverse component of spin, Sg, 
contains the operator e ± ^* e - , which strongly fluctuates 
due to the locking of the dual field </>_ ; therefore the cor- 
relation function (St S^) decays exponentially with the 
distance. Instead, the longitudinal-correlation (SfSf,) 
and the bond ncmatic cor relatiorjDB {Sj Sj +1 Sj, Sj, +X ) 
show power-law decays. Ignoring fluctuations of <p_, we 
obtain the bosonized expressions for these operators as 



Si = M 



47T 



Ai cos 



y^4>+ +tt[M 



sgn(M) 



I- 3 - 
2 



(77) 



+ 1 u2 e WA-R0 + 



2B B 1 e W4 ^ e + cos 



7T / 1 
2 V 2 



\M\ 



(78) 



from which the correlation functions are calculated as 
K + 



{SI St,) = M — 



2tt 2 \£'-1\< 



cos[7r(|Af| - l)(f -£)] 
\V -£\ K +/ 2 



(79) 



(S e S e+1 S e ,S e , +1 ) = B — - - 2/R ^ 



B 



y/ c os[7r(|M| + l)(^ -£)] 

\£l _ £\2/K++K+/2 



(80) 



(73) with B oc A\, B' 



and B" oc B 2 B\. We note 



that the TLL phases with similar power-law correla- 
tions, called the nematic and SDW2 phases, have also 
been discussed for Lhe-roodel ([!]) in a magnetic field-, for 
both ferromagneticEilEjO and antifcrromagnetidlao J\ . 
For small \Ji\/Jq, and A — 1, K + is close to unity, and 
the longitudinal (spin-density-wave; SDW) correlation 
decays more slowly than the nematic correlation. The 
TLL phase with a dominant SDW correlation and short- 
ranged transverse spin correlation is called the SDW2 
state in Refs. fn] and (79[ It is natural to assume that 
the partially polarized phase at A > 1 in zero magnetic 
field is continuously connected to the SDW2 phase in a 
finite magnetic field.Ot§L3 With inter-chain couplings, 
the dominant quasi-long-range SDW correlation, is ex- 
pected to evolve into a true long-rangc-ordcrEJ Since 
K + changes continuously in the TLL phases, it is also 
possible that the system crosses over to a region with the 
dominant nematic correlation (K+ > 2). It is knowri.that. 
such a region does appear at high magnetic fields .OE§E1j 



VI. CONCLUSIONS 

In this paper, we have studied the ground-state prop- 
erties of the one-dimensional spin- g frustrated ferromag- 
netic XXZ model ([!]). In the isotropic case A = 1, the 
nonmagnetic phase in the region —4 < J1/J2 < was 
characterized as the Haldane dimer phase, in which the 
ground state has spontaneous ferromagnetic dimcrization 
and nonlocal string order. We argued that the dimer 
order is associated with an emergent spin-1 degree of 
freedom on every other bond. In the easy-plane case 
< A < 1, the model displays a rich phase diagram 
as m Fig. |. Our previous works have revealed the ap- 
pearance of thergaplcss chiral phase in a wide region for 
—4 < J1/J2 < 0l3 and the unusual alternate appearance 
of the Neel and even-parity dimer phases.cfl In this paper, 
we have newly discovered narrow intermediate gapped 
phases in which the vector chiral order coexists with the 
dimer or Neel order. We described how the properties of 
the various phases can be captured for \Ji\/J2 <C 1 and 
general anisotropy A > by the Abclian bosonization 
formalism, as summarized in Table | (by continuity, the 
same qualitative description can be extended to larger 
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IJilM). 

The Haldane dimer phase we found for A = 1 has 
only a very small excitation gap and, with a weak easy- 
plane anisotropy, is easily replaced by the gapless chi- 
ral phase. With small inter-chain couplings, the gap- 
less chiral phase would evolve into a genuine spiral long- 
range-order. Therefore, the stable appearance of the gap- 
less chiral phase up to the close vicinity of the isotropic 
case A = 1 naturally explains why many quasi-one- 
dimensional cuprates with ferromagnetic J\ < shcpithe 
spiral magnetism and the associated multiferroicity.Ej By 
contrast, it is also expected that the small excitation gap 
(< O.O6J2; see Sec. 1IIB2) in the Haldane dimer phase 
can be enhanced by a coupling with phonons, due to the 
spin-Peierls mechanism as is known inr-the antiferromag- 
netic J1-J2 chain compound CuGeOs.Eil It will be inter- 
esting to explore a spin-Peierls transition to the Haldane 
dimer phase in quasi-lD edge-sharing cuprates without a 
spiral magnetic order. The present study also raises the 
possibility of observing the chiral Haldane dimer state, 
which shows no magnetic order but a spontaneous elec- 
tric polarization due to a vector chiral order of spins. 
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Appendix A: Derivation of the renormalization 
group equations ( |l6| ) 

Here we briefly explain how the RG equations. ( |l6| ) arc 
derived by using the perturbative RG methodcj and the 
operator product expansions (OPE) in the SU(2)i WZW 
theory. 

We first discuss the OPEs in the decoupled spin chains, 
each described by the SU(2)i WZW theory. We drop the 
chain subscript n = 1,2. The OPEs of the uniform spin 
componE nts M il'l obey the well-known SU(2) current 
algcbraHBErH 



M% /L (x,T)M b R/L (0) 



5 ab 



ie abc M R/L {V) 



S7F Z R/L 



2tt; 



R/L 



(Al) 



with zr/l = vr^fix. Here, e abc is the fully antisymmetric 
tensor with e 123 = 1, and summation over repeated in- 
dices are assumed throughout the appendix. The OPEs 
present the singular terms that appear when two opera- 
tors at the points (x, r) and = (0, 0) are brought close 
together. 

The OPEs of the uniform components M R ^ L with the 
staggerecLcomponents N a and the dimerization e are 
given byea 



M a R/L {x,r)N\0) 



M a R/L (x,r)e(0) 



TtW°(0) 



[e o&c A c (0) ± S ab e(0)] 



Attz 



R/L 



(A2) 
(A3) 



These equations imply that M R / L induce mixing of TV 
and e. 

Similar to Eqs. Q^aud (A3), the OPEs among TV 
and e can be derivecEj'CJ by taking advantage of the 
well-known spin-charge separation in ID spin-1/2 Dirac 
fermions; with bosonization, the charge and spin sectors 
of Dirac fermions are described by a free scalar boson 
and the SU(2)i WZW theory, respectively. The use of 
fermionic fields simplifies the calculations of OPEs in the 
WZW theory. For illustration, here we derive the OPE 
of two e's. We take the same conventions as used in the 
Appendix of Rcf. W9l and introduce the right- and left- 



moving fermionic fields ^ r/l,. 
OPEs 



(s =t, I), which obey the 



R/L. 



R/L, 



,(0) 



s ss 



(A4) 



^ Z R/L 

We define the fermionic staggered dimerization operator 



as 



( F = ^ Rs *Ls-^ L ^Rs)- 



(A5) 



Using bosonization, one can show that is related to e 

as 



ep = e cos(v / 2~7r0p), 



(A6) 



where 4> p is the bosonic field of the charge sector. We 
now assume that the charge sector is in the gapped Mott 
phase where 4> p is locked (((f) p ) = 0) as in the Hubbard 
chain at half-filling. This allows us to identify e F with 
Ae, where A = (cos(\/27n/>p)) is a dimensionless constant 
of order unity. The OPE of two e's is then obtained from 
the OPE of two e F 's. 

Performing all possible contractions of four fermion 
fields (see Appendix A of Ref. |82|), the OPE of two ep's 
is calculated as 

e F (x,r)e F (0) = ^ Rs (x,t)V Ls (x,t)¥ Ls ,(0)V Rs ,(0) 



{R O L) 

1 , 1 (pr(Q) 

4tt V 



Pl(0) 



A-kzrzl 
1 



ZL 



zr 



*k s (0)tf Ls (0)*L'(°)*fl.'(0) 
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with 



Pr/l 



t 

R/L, 



,^r/l,s- The last term is related to 



the backscattering term: 

*L*£**L'*««' = ~ 2M R ■ M l - \prPl, (A7) 

where the uniform components of the fermionic spin den- 
sity arc defined as 

Mr = i*L^ *Rs' , M a L = ly Ls a a ss , * Ls , . (A8) 

After gapping out the charge sector, we can neglect the 
fluctuations of Pr/l- Thus we obtain 



0.05 




e(x,r)e(0) 



4tt 2 A< 



ZRZ L 



A 2 



M a (0)-M L (0). (A9) 



Similar calculations yield 



N a {x,T)N b (0) = 



4ir 2 \ 2 z R z L 



is 



2ir\ 2 
1 



M£(0) , M£(0) 



A 2 



iY a (x,r)e(0) 



2ttA 2 



M»(0) M|(0) 



~L 



-TV 



(A10) 



(AH) 



where in Eq. (A.10) is expressed in terms of 

fermionic fields as 



■)ab 

'NN 



1 



(A12) 



For the current purpose, we only need the trace (in 
the spin direction indices) of this term, which gives the 
backscattering term: 0^ N = Mr ■ Ml — \prPl- 

In the limit of weak interchain coupling | J\ | <C J2 , 
the OPEs of the perturbation operators in Eq. ( pi] ) are 
readily obtained from the OPEs of operators in each de- 
coupled chain described above. Given the OPEs, one caa, 
write down the corresponding one-loop RG equations.ES 
For example, if the OPE of marginal operators O a and 
Ob have the form 



O a O b 



A 



a b 



(2n) 2 z R z L 



(A13) 



where \ c ab are dimensionless constants, then the one-loop 
RG equation for the perturbation g c O c has the contribu- 
tion 



dg c 
dl 



9a9bKb 

Attv 



(A14) 



FIG. 17: (Color online) Plots of r 2 \(Sf Sf +r )\ 2 for fixed A = 
and various values of J1/J2 around the TLL-dimer transition 
point studied in Ref. ^j. A logarithmic scale is used for the 
horizontal axis. The symbol "D" indicates the estimate of 
the transition point within the current analysis (with a preci- 
sion of 0.05), whieli agrees reasonably well with the previous 
accurate estimated J2/J1 « 3.0893. 



Appendix B: TLL-dimer transition 



In Sec. 1VB, we determined the BKT transition points 



between gapless and gapped chiral phases by observing 
the logarithmic correction in the spin correlation func- 
tion (Figs. [14] and 15). Here we test the validity of the 



method with a simpler example. We consider the anti- 
ferromagnetic XY model with Ji, Ji > and A = 0. For 
large J\l 3i{^ 3), the system is in a Tomonaga-Luttinger 
liquid (TLL) phase, in wiiich the transverse spin correla- 
tion function behaves asEll 



/ QX QX \ 
\ l D t'l 



\l' - i\n \l' - l\v+i/v 



(Bl) 



Here Aq and A\ are non-universal constants. The decay 
exponent rj gradually increases as J1/J2 is decreased. At 
77 = 1, a BKT transition from the TLL to the singlet 
dimer phase occurs. At the transition point, a multiplica- 
tive logarithmin.correction appears in the spin correlation 
functionSBBEl 



/ QX QX 



\£'-£\ 



■ln5(|£' -£\/a)- 



(B2) 



In Fig. [Tt], we plot the function r 2 \(SfSf +r }\ 2 for various 
J1/J2 around the BKT transition point. From the linear 
behavior as a function of lnr, we locate the BKT transi- 
tion point. In this figure, the data points of J\/ J2 = 3.10 
and 3.15 exhibit almost linear behavior. It is not easy 
to decide which one of the two curves is closer to the 
perfect linear dependence. Here we choose the one with 
smaller correlations since the iTEBD method tends to 
underestimate correlations at large r. The determined 



point J1/J2 = 3.10 agrees reasonably well with the pre- vious accurate estimates J%/ J\ ~ 3.0893 
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